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ABSTRACT 

The origin of interstellar dust in galaxies is poorly understood, particularly the relative con- 
tributions from supernovae and the cool stellar winds of low-intermediate mass stars. Re- 
cently, large masses of newly-formed dust have been discovered in the ejecta of core-collapse 
supernovae. Here, we present Herschel PACS and SPIRE photometry at 70-500 /im of the 
historical, young supernova remnants: Kepler and Tycho; both thought to be the remnants 
of Type la explosion events. We detect a warm dust component in Kepler's remnant with 
Td = 82t^ K and mass ^ (3.1^0 e) ^ 10 '^M©; this is spatially coincident with thermal 
X-ray emission and optical knots and filaments, consistent with the warm dust originating in 
the circumstellar material swept up by the primary blast wave of the remnant. Similarly for 
Tycho 's remnant, we detect warm dust at 90t'^ K with mass (8.6^1 X 10-3 Mq. Comparing 
the spatial distribution of the warm dust with X-rays from the ejecta and swept-up medium, 
and Ha emission arising from the post-shock edge, we show that the warm dust is swept up 
interstellar material. We find no evidence of a cool (25-50 K) component of dust with mass 
> 0.07 M© as observed in core-collapse remnants of massive stars. Neither the warm or cold 
dust components detected here are spatially coincident with supernova ejecta material. We 
compare the lack of observed supernova dust with a theoretical model of dust formation in 
Type la remnants which predicts dust masses of 88(17) x 10"^ M© for ejecta expanding into 
ambient surrounding densities of 1(5) cm^'^. The model predicts that silicon- and carbon- 
rich dust grains will encounter, at most, the interior edge of the observed dust emission at 
~400 years confirming that the majority of the warm dust originates from swept up circum- 
stellar or interstellar grains (for Kepler and Tycho respectively). The lack of cold dust grains 
in the ejecta suggests that Type la remnants do not produce substantial quantities of iron-rich 
dust grains and has important consequences for the 'missing' iron mass observed in ejecta. 
Finally, although, we cannot completely rule out a small mass of freshly-formed supernova 
dust, the Herschel observations confirm that significantly less dust forms in the ejecta of Type 
la supernovae than in the remnants of core-collapse explosions. 

Key words: Supernovae: individual: Kepler, Tycho - ISM: submillimetre dust: ISM - Galax- 
ies: abundances - submillimetre 
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1 INTRODUCTION 

Dust in galaxies is thought to be produced by in the stellar winds 
of both low-intermediate mass (LIM) Asymptotic Giant Branch 
(AGB) stars (e.g. Gehrz 1989; Whittet 2003; Sargent et al. 2010) 
and, to an unknown extent, by massive stars (Hofner 2009; Gomez 
et al. 2010; Andrews et al. 201 1 ; Gall, Hjorth & Anderson 2011a). 
Massive stars may also produce dust in the expanding gas ejecta 
when they explode as supernovae (e.g. Clayton et al. 2001; Todini 
& Ferrara 2001). In the Milky Way (see Zhukovska & Gail 2008), 
the major dust source is presumed to be LIM stars, but when ac- 
counting for dust destruction timescales and the observed dust mass 
in the Galaxy, the dust injection rate from stars required is an order 
of magnitude higher than observed. An alternative source of dust is 
required to make up the dust budget (e.g. Jones 2001). This shortfall 
in the dust mass estimated from the dust injection rates from LIM 
stars is also observed in the interstellar medium (ISM) of the Large 
Magellanic Cloud (Matsuura et al. 2009), in galaxies out to z = 0.4 
(Dunne et al. 2011) and in dusty high-redshift galaxies (see Mor- 
gan & Edmunds 2003; Dwek, Galliano & Jones 2007). Although 
recent work by Jones & Nuth (2011) suggest that the dust-budget 
problem is not as problematic as previous works suggest given the 
uncertainties in mass loss rates and the efficiencies of dust forma- 
tion in stellar winds, it is clear that significant dust production in 
SN ejecta would alleviate this dust budgetary problem. The aver- 
age dust yield per SN required to explain dusty galaxies at both low 
and high redshifts is of order 0.5 — 1 M© (Michalowski, Watson & 
Hjorth 2010; Gall, Anderson & Hjorth 201 lb). 

Although SNe have long been proposed as a source of dust in 
the ISM (Dwek & Scalo 1980; Clayton et al. 2001), observational 
evidence has been scarce until recently. Wooden et al. (1993) first 
detected emission from 10"* Mq of dust in SN1987A using the 
Kuiper Airborne Observatory. Subsequent Mid-far-infrared (FIR) 
observations of Galactic and extra-galactic core-collapse remnants 
with the Spitzer Space Telescope detect small quantities of warm 
dust(~ (1-50) X 10"^ M©) in young (t < 1000 days - Sugerman 
et al. 2006; Meikle et al. 2007, 201 1; Kotak et al. 2009; Andrews et 
al. 2010, 2011; Fabbri et al 2011; Szalai et al. 2011) and old rem- 
nants (Williams, Chu & Gruendl 2006; Rho et al. 2008). The evi- 
dence for dust production in Type la remnants is scarce (Borkowski 
et al. 2006; Seok et al. 2008): mid-IR observations of the Kepler 
SNR with the Infrared Space Observatory (ISO; Douvion, Lagage 
& Pantin 2001) found thermal emission from warm dust (100 K) 
with mass 10~* M©. More recently, Ishihara et al. (2010) detected 
10"^ M0 of swept-up (warm) dust in the Tycho SNR with AKARI, 
with a tentative suggestion that 10~* M0 of dust seen in the north- 
east of the remnant may have been formed in the ejecta. Outside 
of our Galaxy, Borkowski et al. (2006) obtained mid-IR observa- 
tions of four la remnants in the Large Magellanic Cloud, finding 
< 10~^ M0 of swept-up interstellar dust. The dust masses from 
the mid-far IR studies are orders of magnitude lower than needed 
to solve the dust budget issues. 

In order to investigate if supernova remnants (SNRs) contain 
cold dust missed by previous infrared (IR) observations, the Sub- 
millimetre Common User Bolometer Array (SCUBA) was used 
to observe Cassiopeia A (Cas A), a remnant of a Galactic core- 
collapse SNe (Krause et al. 2008a). The cold dust seen at 450 and 
850 /im was interpreted as emission from dust associated with the 
remnant due to the spatial correlation with the X-ray, radio and sub- 
millimetre (submm) emission (Dunne et al. 2003). Subsequently, a 
significant fraction of the submm emission was shown to originate 
from molecular clouds along the line of sight (Krause et al. 2004; 



Wilson & Batrla 2005). Recent observations with the Balloon- 
borne Large Aperture Submillimeter Telescope, BLAST and the 
Herschel Space Observatory showed it was difficult to distinguish 
between cold dust in the remnant and cold dust from intervening 
interstellar clouds using photometry alone (Sibthorpe et al. 2010; 
Barlow et al. 2010). These studies revealed a new cool component 
of dust in Cas A with mass 0.08 Mq at ^ 35 K, yet even if all 
of the dust survived the passage through the shock, the dust mass 
is about one order of magnitude lower than necessary to solve the 
dust budget problem. 

A possible method to distinguish between supernova dust and 
contamination from interstellar material along the line of sight was 
proposed by Dunne et al. (2009) where they traced the alignment 
of the dust towards Cas A with the SNR magnetic field using po- 
larimetry at 850 ^m. They find that a significant fraction of the 
submm emission towards Cas A is strongly aligned with the mag- 
netic field of the remnant. The dust mass inferred from the po- 
larimetry study of Cas A agrees with the dust masses estimated 
from Herschel photometry observations of the core-collapse SNe 
SN1987A (Matsuura et al. 2011); in the former 0.2 - 1.0 M© of 
dust was detected depending on the composition. The limit on the 
ejected heavy element mass in these sources places a further con- 
straint on the maximum physical dust mass, finding 0.2 < Md < 
0.5 Mq. 

The evidence for dust formation in core-collapse SNRs is ac- 
cumulating, though it is not yet understood whether the amount of 
SN dust created is dependent on the progenitor mass, the metallic- 
ity or the type of explosion. Neither is it clear how much dust will 
survive, though theoretical studies suggest a significant fraction of 
dust grains formed in ejecta will be destroyed before injection into 
the ISM (Jones 2001; Bianchi Sl Schneider 2007) particularly in 
the ejecta of envelope-poor core-collapse SNe or Type la remnants 
(Nozawa et al. 2010; Nozawa et al. 201 1). As yet, there is no clear 
evidence for freshly formed dust in the remnants of Type la explo- 
sions, though theoretical models of dust formation in the remnants 
of Type la events suggests that up to 0.2 Mq dust could form in the 
ejecta (Nozawa et al. 201 1). 

Historical remnants are unique probes for understanding dust 
production and indeed, destruction following core collapse Type II 
supernovae as well as low mass binary Type I events. In the Her- 
schel era, we are finally beginning to probe the contribution of SNe 
to the dust budget, yet given the somewhat large uncertainties in our 
understanding of the dust yield from stellar sources in other galax- 
ies and in interpreting statistical descriptions of the dust content in 
galaxies from wide-area surveys (e.g. Dunne et al. 2011), it is im- 
portant to investigate dust production in those nearest, youngest and 
resolved galactic remnants. One can also argue that understanding 
dust formation in supernovae is important regardless of the explo- 
sion mechanism: dust formation following Type-II SNe might sug- 
gest that Type-la SNe would also be likely dust producers (Clayton 
et al. 1997; Travaglio et al. 1999). Alternatively, detecting dust from 
circumstellar material swept up by the blast wave of a Type la can 
be a useful probe of the progenitor and explosion mechanism. 

The resolution, wavelength coverage and sensitivity of Her- 
schel allows us to investigate the origin and quantity of dust in 
historical supernova remnants. Here we present IR and submm ob- 
servations (iQ of the Tycho and Kepler remnants. We determine 
the dust mass in and towards the remnants in ^and ^and com- 
pare the IR emission with multiwavelength tracers in O.ll and i]4. II 
The origin of the dust is discussed in ^and we use the theoretical 
model of Nozawa et al. (201 1) to compare the expected dust masses 
with the observations. The results are summarised in ^ 
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2 OBSERVATIONS AND DATA REDUCTION 

2.1 Herschel Photometry 

The Kepler and Tycho supernova remnants were observed with the 
Herschel (Pilbratt et al. 2010) Photodetector Array Camera and 
Spectrometer (PACS; Poglitsch et al. 2010) and Spectral and Pho- 
tometric Imaging Receiver (SPIRE; Griffin et al. 2010) at 70, 100, 
160, 250, 350 and 500 /im as part of the Mass Loss from Evolved 
StarS survey (MESS: Groenewegen et al. 2011). Kepler was ob- 
served on 10th Sept and 1 1th March 2010; Tycho on 20th March 
2010 and 18th Jan 2010 for PACS and SPIRE respectively. The 
PACS data were obtained in 'scan-map' mode with speed 20"/s in- 
cluding a pair of orthogonal cross-scans over an area ~ 22' x 22'. 
The SPIRE maps were observed in 'Large Map' mode with scan 
length of 30' over ~ 32' x 32'; a cross-scan length is also taken, 
with a repetition factor of three. The data were processed following 
the detailed description given in Groenewegen et al. (201 1). 

The PACS maps were reduced with the Herschel Interactive 
Processing Environment (HIPE; Ott 2010) applying all low-level re- 
duction steps (including deglitching) to Level 1 . The SCANAMOR- 
PHOS software (Roussel 2011) was then used remove effects due 
to the thermal drifts and the uncorrelated 1 // noise of the individ- 
ual bolometers. Final Level 2 map creation through data projection 
onto a changeable spatial grid was performed with SCANAMOR- 
PHOS. Since the HIPE reduction was based on the calibration file 
version v5, interim calibration factors of 1.12, 1.15 and 1.17 were 
applied to the PACS data at the respective wavelengths to be consis- 
tent with the current calibration file v6. The Full Width Half Max- 
imum (FWHM) at 70, 100 and 160 /im is 6, 8 and 12" respectively. 
The flux calibration uncertainty for PACS is currently estimated as 
10 per cent for the 70 and 100 /im bands and 20 per cent at 160 ^m 
(Poglitsch et al. 2010). Colour corrections for PACS are expected 
to be small compared to the calibration errors. To correct for colour 
with dust spectrum v^'^ and temperatures of 100 K, we divided by 
1.00 and l.o£]at 70 and 100 ^m. The 20 K interstellar dust compo- 
nent dominates the 160 /im image, and as such the correction factor 
is 0.957. 

For SPIRE, the standard photometer pipeline (HIPE v.5.0) was 
used (Griffin et al. 2010) with an additional iterative baseline re- 
moval step (e.g. Bendo et al. 2010). The SPIRE maps were cre- 
ated with the standard pipeline NAIVE mapper. We multiply the 
350 /^m data product by 1.0067 to be in line with most recent cali- 
bration pipeline (equivalent to HIPE v7 reduction). Since the beam 
areas increase as a function of pixel scalfl the FWHM is 18.1, 
24.9 and 36.4" for pixel sizes of 6, 10 and 14", at 250, 350 and 
500 respectively. The SPIRE calibration methods and accura- 
cies are outlined by Swinyard et al. (2010) and are estimated to 
be approx 10 per cent. The pipeline produces monochromatic flux 
densities for point sources with vSu oc u^'" but at the longer 
wavelengths, colour corrections become significant. We multiplied 
by the colour corrections 0.999, 1.0045, 1.0302 as defined in the 
Herschel SPIRE manual^, appropriate for /3 = 1.5 and extended 
sources. 

2.2 Multiwavelength Datasets 

Optical data for the Tycho and Kepler remnants were obtained 
from the Hubble Space Telescope (HST) archive via ALADIN 

^ http://herschel.esac.esa.int/twiki/pub/PubUc/PacsCaUbration'Web/\ 
^ http://herschel. esac. esq. int/Docs/SPIRE \ 



(Bonnarel et al. 2000). For Kepler's remnant we used observations 
in the filters 658, 660 and 502 nm where the emission originates 
from Hof + [Nil], [Nil] and [Olll] respectively (e.g. Sankrit et al. 
2008). Bad pixels were masked by hand. For Tycho, we use data 
observed at 656 nm (Lee et al. 2010), tracing Ha in the eastern 
limb of the remnant. 

Spitzer Multiband Imaging Photometer (MIPS; Rieke et al. 
2004) Level 2 calibrated data were obtained via the archiv^f] (Ke- 
pler: PI W.Blair; Tycho: PI J.Rho) at 24 and 70 ^m. For 100 ^im 
data, we use the Improved Reprocessing of the IRAS Survey (IRIS; 
Milville-Deschenes &. Lagache 2005) which have been calibrated 
with respect to data taken with the Diffuse Infrared Background 
Experiment (DIRBE). Wide field Infrared Survey Explorer (WISE; 
Wright et al. 2010) data were also obtained at 22 pm. 

The radio data at 18 cm for Tycho was obtained through the 
Very Large Array (VLA) archive. The Kepler data at 6 cm was 
kindly provided by T. DeLaney (DeLaney et al. 2002). The Chan- 
dra X-ray data for Kepler (PI S. Reynolds; Reynolds et al. 2008) 
and Tycho (PI J.Hughes; see Badenes et al. 2006; Katsuda et al. 
2010) were available online as calibrated files extracted into soft 
and hard regions of the X-ray spectrunjfl For Kepler's remnant, the 
X-ray data used here range from 0.30 - 0.72, 0.72 - 0.90 and 
1.7 - 2.1 keV, and for Tycho, we have 0.95 - 1.26, 1.63 - 2.26 
and 4.1 — 6.1 keV. These ranges show line emission from silicon 
and iron orginating from the ejecta as well as soft and hard (with 
thermal and non-thermal contributions) X-ray continuum orginat- 
ing from the blast wave. 

To compare the submm data with the molecular emission 
towards and surrounding the remnants, we use high resolution 
^^CO(J = 2-1) and "C0( J = 2-1) observations for Ke- 
pler (see Gomez et al. 2009 for full details). We also include our 
SCUBA image of Kepler at 850 ^m (Morgan et al. 2003; Gomez et 
al. 2009). Low resolution carbon monoxide '^^C0( J = 1 — 0)maps 
for Tycho's remnant were taken from the Canadian Galactic Plane 
Survey (CGPS; Taylor et al. 2003). The data were integrated over 
—68 to — 53kms^^ to compare the remnant with molecular gas 
thought to be interacting with the supernova remnant (see ^STTJ. 

2.3 Zero Levels 

Due to the difficulty in finding a suitable baseline with which to 
subtract the background in the Herschel SPIRE maps, we correlate 
the pixel intensity from Herschel with the archival IRAS IRIS maps 
(calibrated with respect to DIRBE) to estimate the zero level (see 
Fig. [T] Table [T). The Herschel data were convolved to the same 
angular resolution as the IRIS 100-/im map (258") and regridded 
to 90" pixels. A least squares fit analysis was performed on the 
Herschel vs. IRIS data (whilst masking out the remnant) to de- 
termine the flux offsets; this sets the zero level of the PACS and 
SPIRE data assuming that the IRIS data is a true representation of 
the sky. A further scaling factor needs to be accounted for due to 
the expected colour correction between the bands when compar- 
ing to the 100 pva band (the so-called 'gain' factor). These were 
estimated by interpolating a modified blackbody curve u^B{iy, Td) 
(after a preliminary correction of the zero levels). The zero level 
of the PACS and SPIRE maps are negative except for the PACS 
70 pm map of Kepler. The least squares fit was then repeated (see 
e.g. Bracco et al. 2011). The errors on the fitting were estimated 

^ http://irsa.ipac.caltech.edu/data/SPITZER/docs/spitzerdataarchives/\ 
^ http://chandra.harvard.edu/pwarmo/openFITS/\ 
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Wavelength Zero level (MJy/sr) 

{p,m) Kepler Tycho 



70 


4.7 ±0.4 


-0.9 ±0.1 


100 


-3.1 ±0.2 


-2.5 ±0.3 


160 


-4.8 ± 0.3 


-2.6 ±0.2 


250 


-37.4 ±3.8 


-30.8 ±4.8 


350 


-30.5 ±2.0 


-20.4 ±2.1 


500 


-25.2 ±1.4 


-7.7 ± 1.3 



50 p 
40 - 
30 - 




IRIS 100 (MJy/sr) 

Figure 1. The Herschel SPIRE intensities for the Kepler observations ver- 
sus IRAS IRIS 100-/im in MJy/sr. The SPIRE maps are convolved to the 
same resolution as the IRIS maps and regridded onto the same pixel scale. 
The red line is the least squares function, the grey circles are the Herschel 
data divided by the gain factor (due to the colour correction between IRIS 
at 100 /.tm and the SPIRE wavelengths). The data are masked at the location 
of the remnant. The grey shaded regions indicate the results from the 1000 
bootstrap fits to the same dataset. 

using the confidence intervals from the bootstrap method, where 
1000 fits were made to the bootstapped data, Table [T] This error is 
added in quadrature along with the calibration errors listed above 
and introduces (at most) an error of ^ 17% (well within the 20% 
calibration error used for PACS at 160 /im); note that this does not 
account for any systematic errors. 



3 KEPLER'S SUPERNOVA REMNANT 

The Kepler supernova event was first seen by Johannes Kepler in 
1604 and the remnant has a shell-like structure ~3 arcmin in diam- 
eter. Estimates of its distance, using HI absorption features, range 
from 3.9 to 6kpc (Reynoso & Goss 1999; Sankrit et al. 2008). Its 
classification has been controversial (see the excellent reviews in 
Blair et al. 2007; Reynolds et al. 2008; Sanki-it et al. 2008) with 
evidence pointing towards both Type-la - the thermonuclear explo- 
sion of a low-mass accreting star in a binary system - or Type-lb - 
the core collapse of a massive star. The deepest X-ray analysis with 
Chandra (Reynolds et al. 2008), however, suggests that a Type-la is 



Table 1. Zero level on the PACS and SPIRE maps calculated using least 
squares fit to Herschel data and IRAS IRIS data. Errors are the 1 cr estimated 
from bootstrapping the data 1000 times and refitting each set. 




Figure 2. Three-colour FIR-submm image of Kepler's SNR with PACS at 
70 (blue), 100 (green) and 160 /im (red). The image has been smoothed. 
The area shown is 16.5' X 16.6'. 



the most likely classification using both iron abundance constraints 
and the lack of a neutron star. The deep X-ray images also provide 
evidence of the blast-wave interacting with dense circumstellar ma- 
terial (CSM). This is unexpected for Type la events although ob- 
servations of CSM are becoming more widespread; the presence of 
such a dense CSM surrounding Kepler led Reynolds to suggest the 
thermonuclear explosion of a single massive star (Af« ~ 8Mq) 
as an alternative to the canonical White-Dwarf (WD) binary sce- 
nario of la explosions. Whether Kepler was the result of a single 
progenitor or WD-AGB-binary system explosion, the mechanism 
was likely to be a deflagration explosion (Yang et al. 2009; see 
Woosley, Taam & Weaver 1986 for a review on explosion models). 
In this scenario, the heating of the ejecta by radioactive elements 
is higher than in core-collapse remnants, the ejecta densities will 
also decrease rapidly, though the ejecta is well mixed since the re- 
verse shock progresses quickly through the layers in a deflagration 
explosion (see H5.l\ . 

The warm dust emission detected in Kepler with ISO (Dou- 
vion et al. 2001), follows the well-known north-south shell asym- 
metry seen also at optical. X-ray and radio wavelengths. The cor- 
relation of the dust with the optical Ha suggested the dust was 
shock-heated circumstellar dust, swept up by the Shockwave of the 
blast. Blair et al. (2007) used Spitzer to revise the dust mass to 
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5.4 X 1O~*M0. They also found that emission at 70 /xm arises 
only in the brightest knots of the shell seen at 24 /xm. The 160 /im 
Spitzer data was poor-quality due to its low resolution and strong 
background gradient. Morgan et al. (2003) originally interpreted 
their longer wavelength 450 and 850-/im SCUBA data as evidence 
for 0.3-1 M© of cold dust associated with the remnant. They com- 
pared the emission from interstellar material towards Kepler us- 
ing submm continuum imaging and spectroscopic observations of 
atomic and molecular gas, via HI, ^^CO and "CO and detected 
weak CO emission from diffuse, optically thin gas at the locations 
of some of the submm clumps. The contribution to the submm 
emission from foreground molecular and atomic clouds using these 
CO tracers was found to be negligible, but these emission lines may 
not be tracing all of the gas in the ISM. 

The three colour Herschel PACS image of Kepler is displayed 
in Fig.|2] The remnant is clearly seen in blue and follows the same 
asymmetric structure as seen in the X-ray and radio (i ]3.1t : the 
emission is brightest in the north along an arc-like shell and some 
emission is seen towards the centre of the remnant. The low level 
red structures seen across the map (outside of the remnant) origi- 
nate from dusty interstellar clouds with temperatures ~ 13 — 20 K. 



3.1 The multiwavelength view 

In Fig. [3] we present a multiwavelength view centred on Kepler's 
SNR, comparing the Herschel PACS and SPIRE photometric data 
with synchrotron emission from the remnant at 6 cm, unrelated 
molecular gas (traced by the CO emission) and 850 emission 
from cold dust. The complete shell-like structure of SN material is 
clearly seen at radio wavelengths, tracing the forward and reverse 
shocks (DeLaney et al. 2002). In Fig. |4]we compare the optical 
and X-ray emission with the warm dust from Herschel-PACS at 
70 /xm. The well-known north-south asymmetry seen at optical, ra- 
dio and X-ray wavelengths with the enhanced northern emission 
(thought to originate from the motion of Kepler through the sur- 
rounding interstellar material) is also seen in the PACS data. The 
'arm' structure across the middle (clearly seen in the radio. Fig. [3} 
originates from a projection effect of material in the front and back 
of the shell (Blair et al. 2007). Herschel detects faint IR emission 
arising from this component at 70 /im. The faint 'ear' -like features 
seen along the outermost midsection of the expanding shell at ra- 
dio and X-ray wavelengths (Figs. [3] & are not seen in the in- 
frared. At the longer (160 nm) PACS and the SPIRE wavelengths, 
it is difficult to disentangle any SNR emission from the large-scale 
interstellar structures extending across the region; the maps are in- 
creasingly dominated by unrelated dust clouds at temperatures of 
~ 16 — 20 K. The southern cloud of dust which is outside the area 
of the blast wave and ejecta regions is the same cloud orginally de- 
tected in molecular emission by Gomez et al. (2009) at 5 kms~^. 
The cold dust cloud in the north-east is coincident with the CO 
structure at 10 km (see i |5.3| for more discussion). 

Optical emission arising from dense, knotty regions in the cir- 
cumstellar medium overrun by the blast wave, is traced in [Nil] 
and Hq. Interestingly, the 70 fim emission correlates with the Ha 
-I- [Nil] features (Fig.|4]l in both the northern limb (particularly in 
the north-west region) and across the projected shell along the mid- 
dle. Since the Ha emission originates from coUisional heating of 
the thin region behind the SN shockfront (i.e. from swept-up sur- 
rounding material), the spatial agreement with the PACS emission 
and the ionised features suggest this dust component arises from 
the swept up medium. This is further supported when comparing 
the FIR with X-ray emission. There are three X-ray components 



in Kepler: continuum X-rays from the blast wave sweeping up sur- 
rounding circumstellar and interstellar material, X-ray line emis- 
sion from the warm ejecta, and continuum emission arising from 
the non-thermal synchrotron emission in the north, south and east- 
em limbs. The warm dust emission in the FIR is confined within 
the radius of the outer blast-wave and is seen only at the brightest 
points in the optical, radio and X-ray (Fig. |4] - as noted in Blair 
et al. 2007) with peaks at 70 /im seen where the soft and hard ther- 
mal X-rays overlap. Therefore, the warm dust is associated with the 
densest gas behind the forward shock rather than the inner ejecta 
material; this confirms the dust originates from the outermost re- 
gions of the remnant i.e. from the swept up shocked circumstellar 
and/or interstellar medium ( i]5.2t . 

3.2 The spectral energy distribution 

The IR-submm spectral energy distribution (SED) towards Kepler's 
SNR is shown in Fig. [5] The fluxes estimated from Spitzer, Her- 
schel and SCUBA in this work within an aperture of radius 120" 
centred on the remnant, are plotted with their errors (estimated from 
the calibration error, the error in the zero level estimated in i ]2.3| and 
from the variations in the pixel RMS within the aperture). All the 
data were smoothed to the resolution of the SPIRE 500 pm data 
using a Gaussian kernel. The aperture is somewhat larger than the 
SNR boundary defined by the edge of the radio (103") but was 
chosen to encompass the SNR emission in the images smoothed 
to the 500 fim beam. The Spitzer fluxes were colour corrected us- 
ing algorithms detailed in the MIPS handboolB We also compare 
our fluxes with previously published fluxes in the literature (Arendt 
1989, Braun 1987, Douvion et al. 2001, Morgan et al. 2003, Blair et 
al. 2007) although these were obtained using apertures with radius 
ranging from 100 — 120" and therefore differ from the fluxes in Ta- 
ble|2]accordingly. Note also that we measure the flux encompassed 
in the aperture around Kepler, and therefore includes emission from 
the foreground cloud. 

The contribution to the IR/submm fluxes from synchrotron 
emission is obtained by scaling the VLA radio map of Kepler us- 
ing a constant spectral index. Si, oc i/" where the mean value of 
a is -0.71 but ranges from -0.6—0.85 (DeLaney et al. 2002). The 
synchrotron emission at 500 pm is 0.20 ± 0.06 Jy (~ 15 percent 
of the total flux). We fit a two-component modified blackbody to 
the dust SED after the synchrotron component has been subtracted; 
the model is the sum of two modified Planck functions each with a 
characteristic temperature, r„ and Tc (Eq.[T](: 

S,^N^xu^B{u,T^)+NcXu^B{i^,Tc) (I) 

where A*'„ and Nc represent the relative masses in the warm 
and cold component, B{i/,T) is the Planck function and P is the 
dust emissivity index. The model was fitted to the SED (constrained 
by the 12 pm flux) and the resulting parameters (A^,„, Nc, T^, Tc, 
/3) which gave the minimal chi-squared were found. We also ap- 
plied a bootstrap analysis to our SED-fitting to determine the errors 
on the SED model. The photometry measurements were perturbed 
according to the errors at each wavelength, the new fluxes were fit- 
ted in the same way, and the SED parameters recorded, with the 
procedure repeated 1000 times. 

The total dust mass is calculated using Eq.[2l 



Md = X 

«:500 



+ 



Nc 



B(500 ^m, r„) B(500 ^m, T, 



(2) 



\http://irsa. ipac. caltech. edu/data/SPITZER/docs/mips/\ 
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Figure 3. Muldwavelength montage of Kepler's supernova remnant and its environs. Top: PACS and middle: SPIRE bands centred at RA = 262.671°, 
Dec = —21.4914° (J2000.0; yellow cross); the region show is 7.2' X 7.2'. The location of the forward shock (radius 103") is indicated with the white circle 
(DeLaney et al. 2002). Also shown are ioffom: the SCUBA 850 ^tm (Morgan et al. 2003); theCO(J = 2 - 1) map integrated from -200 <v < 200 km s"! 
(Gomez et al. 2009) and the 6 cm VLA map (DeLaney et al. 2002). 









IR Fluxes Towards Kepler SNR (Jy) 








24/xm 


70/im 


100/im 


160/im 


250^m 350/^m 


500Afm 


850Atm 


Published 


9.5 ± 1.0° 


5.6 ± 1.4'' 


2.9 ±1.1* 


< 0.9"^ 




3.0 ±0.2=-'^ 


1.0 ±0.2'= 


This work 


9.5 ± 1.0"= 


12.3 ±2.7 


11.2 ± 2.4 


16.5 ± 2.9 


13.0 ±2.8 5.8 ±1.2 


2.7 ±0.6 


0.7±0.1 


Nonthermal 


0.02 


0.05 


0.06 


0.09 


0.12 0.16 


0.20 


0.29 



Table 2. Fluxes in Jy measured for Kepler's remnant within the 120" aperture. Note that the previous published integrated fluxes use different aperture sizes 
to this work. References: " - Blair et al. (2007); - average values from Douvion et al. (2001), Arendt (1989) and Saken, Fesen & Shull (1992) IRAS data; " - 
Gomez et al. (2009); - 450 /im. " - Measured on Spitzer MIPS archive data. 



where K500 is the dust absorption coefficient and D is the distance 
(we choose the 500 /xm flux since this is less sensitive to the temper- 
ature). The distance to the remnant is uncertain, with measurements 
of Hi absorption placing Kepler at ~ 5 kpc (see Reynoso & Goss 
1999 for a full discussion). Sankrit et al. (2008) have revised the 
measurement to 3.9!ti'4 kpc using HST proper motion data. We 



follow Blair et al. (2007) in adopting a distance of 4 kpc. To deter- 
mine the dust mass, we use K500 = 0.1 kg~^ with k oc , 
appropriate for grains in the Milky Way and j3 ^ 1 ~ 2 (Draine 
& Lee 1984). The best-fit dust mass within the aperture obtained 
from the SED model (Eqs.E&|2j is therefore M„ = 0.0035 M© 
and Mc = 2.2 M©, with temperatures T„ = 84 K, = 19 K and 
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Figure 4. Compaiison of the warm dust emission traced by PACS 70-/xm data (wliite contours) witli optical and X-ray. Left: Three colour image of Kepler's 
supernova remnant with the HST at 658 nm (red), 660 nm (green) and 502 nm (blue) indicating emission at Ho + [Nil], [Nil] and [OlII] respectively. The 
optical emission mostly arises from shocked circumstellar/interstellar material interacting with the primary blast wave. The north component indicated with 
the fiducial is thought to be ejecta material (Sanrkit et al. 2008) which has overrun the primary shock boundary due to the contact discontinuity in the blast wave. 
Right: Three colour X-ray image using Chandra archive data showing soft (0.3 ~ 0.72 IscV: red), (0.72 — 0.9 IteV: green) and hard X-rays (1.7 — 2.1 keV: 
blue). The images are centered on aj2000 = 262.663° and (5j2ooo = —21.485° and are 3' X 3' across. 




Wavelength C/im) 

Figure 5. The IR-radio SED towards Kepler's SNR. The previously pub- 
lished fluxes from the literature (Table |2) are shown with errors as ver- 
tical bars (Spitzer = purple triangles; IRAS = blue diamonds; SCUBA = 
orange squares). The solid black line is the sum of the two-temperature 
thermal dust components fitted to the IR (T^ = 86 K and Tc = 20 K) 
and submm fluxes (dot-dashed and dashed lines respectively). The dot-dot- 
dot-dashed line shows the non-thermal SED expected from the synchrotron. 
Si, oc jy^O.Ti scaled by the 6 cm fluxes. The red line shows the total SED 
from the non-thermal and thermal components. Note that the previous pub- 
lished integrated fluxes use different aperture sizes to this work. 



P — 1.5. The temperature of the warm component is consistent 
with coUisional heating of dust in the shocked gas with tempera- 
ture and electron density of 1.2keV and ris ~ 20 cm""^ (Bouchet 
etal. 2006, their Fig. 15). 

The median parameters and confidence intervals from the 
bootstrap method are Hsted in Table |4] The dust mass depends 
strongly on the choice of k, for example, typical interstellar dust 
grains are described well with k — 6.7m^lcg~^ at 10 fim (e.g. 
Draine & Li 2001). In this case, the warm dust mass would be re- 
duced to AIui = 1 X 10~^ Mq. We note that the uncertainty in k 
could be an order of magnitude. 

Temperature maps of the dust emission towards Kepler were 
created using the convolved IR-submm maps from 24-500 /im, and 
fitting a two-temperature modified blackbody SED to each pixel. 
We extracted the warm and cold components for each pixel, shown 
in Fig.|6](a) and (b). In order to account for the non-thermal compo- 
nent, the flux due to synchrotron radiation was subtracted from the 
dust maps on a pixel-by-pixel basis before the greybody fits were 
made. The synchrotron contribution in each pixel was estimated 
using a spectral index map created from the 6 cm and 23 cm VLA 
observations convolved to the same beam (see also DeLaney et al. 
2002). The somewhat flatter spectral index in the northern bright 
ring (a ~ —0.68 compared to the average value of -0.71) resulted 
in removing most of the cold dust in this region, though there is 
residual emission associated with this component within the SNR 
shell (Fig.|6](b)). This cold dust feature appears to be coincident 
with a similar structure seen in the warm dust component but the 
cold dust is further out from the centre by ^ 10" compared to the 
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(b) 



Figure 6. Dust temperature maps derived from fitting a modified greybody 
to file 24-500 ^m maps towards Kepler's supernova remnant sliowing(a) 
warm dust and (b) cold dust. The colourbar shows the temperature in K. The 
solid black circle indicates the photometry aperture used with radius 120" 
and the black dashed circle is the average radius of the outer shockfront 
103". 



warm dust emission which is spatially coincident with the radio 
(Fig.[6](a)). 

The cold dust structures in the north are also coincident with 
the structures seen at 850 /^m with the SCUBA camera (Gomez et 
al. 2009, Fig. [Sj originally attributed to supernova dust due to the 
lack of correspondance between the CO gas from molecular inter- 
stellar clouds towards Kepler (Fig. O and the 850 /im emission. 
The larger Herschel map clearly shows large-scale and exended in- 
terstellar dust structures right across the entire region towards Ke- 
pler and implies that that the argument in Gomez et al. was flawed. 
This could simply be due to a non-constant ratio of molecular gas- 
to-dust, or that the '^^CO and ^^CO lines are not effective tracers 
of the gas in this region. Alternatively, it is possible that, towards 



Kepler, the interstellar dust structures are associated with atomic 
rather than molecular gas. 

If the cold dust in the north of the remnant is located at the 
distance of the remnant, the location of the emission is suggestive 
of dust swept up by the blast wave, though the mass of dust in the 
northern shell is rather substantial at ~ 1.0 Mq for T ~ 20 K. 
Elsewhere, the cold dust component does not spatially coincide 
with any of the SNR tracers including X-ray, radio and optical 
emission and is unlikely to be related to the SN itself. We discuss 
the origin of the cold dust further in i]5.3l 

Finally, we also tested for the presence of cool dust with tem- 
peratures in the range 30-40 K as found in Cas A (Sibthorpe et al. 
2010; Barlow et al. 2010): the flux contribution of the warm com- 
ponent (visible at 24 ^m) was subtracted from the 70 and 100 ^m 
maps, however, we find no evidence of a cool dust component. 



4 TYCHO'S SUPERNOVA REMNANT 

The Tycho supernova event was observed by Tycho Brahe in 1572. 
It lies at a distance of 1.5 — 3.8 kpc (Reynoso et al. 1997; Krause 
et al. 2008b) and has a thin shell-like morphology with diameter 
~ 8' across. Here, we use the distance given in D = 3.8 kpc as 
this is consistent with the historical records of the explosion, the 
non-detection of the remnant and the distance to the proposed bi- 
nary companion (Lu et al. 2011) (see Krause et al. 2008b for more 
details). The outer blast wave extends to an average distance of 
251" from the centre, with the contact discontinuity between the 
reverse shock material and swept up ISM at angular distances of 
183, 241 and 260" depending on the position angle (Warren et al. 
2005). Tycho's light echo spectrum has revealed it to have been a 
Type la event (Krause et al. 2008b); it is the remnant of a canoni- 
cal single-degenerate binary explosion (Lu et al. 2011) as evident 
from X-ray ejecta abundances (Decourchelle et al. 2001). This is 
further supported by the lack of similarity in Tycho's optical light 
curve with typical sub-luminous (van den Bergh 1993) and over- 
luminous (Ruiz-Lapuente 2004) Type la explosions (Krause et al. 
2008b). Badenes et al. (2006) performed detailed comparison be- 
tween the X-ray spectra and the results from hydrodynamical mod- 
eling of the supernova shocks, and found that the delayed detona- 
tion explosion model was the best fit to the data. In this scenario, 
the SNR is expanding into an ambient density of ~ 0.6 — 3 cm~^, 
with an ejecta mass of 1.3 M©. 

Many authors have suggested that Tycho is interacting with 
molecular clouds (Reynoso & Goss 1999; Lee, Koo & Tatenatsu 
2004; Cai, Yang & Lu 2009; Xu, Wang & Miller 2011) partic- 
ularly with a large CO structure seen in the north and north-east 
at velocities between -63.5 and -61.5 kms^^; this structure is be- 
lieved to be responsible for the asymmetry seen in the X-ray and 
radio (see Fig. [8]l. Cai et al. (2009) showed that interaction with 
the ISM occurs along the entire boundary of the remnant with in- 
terstellar structures detected at velocities -69 — -58kms~^; the 
total gas mass of clouds interacting with Tycho is 10'^ M0 (Xu et 
al. 2010). Conversley, Tian & Leahy (2011) propose that the CO 
clouds are not interacting with the remnant since the densities of 
the molecular gas clouds in the north are ^ 200 cm"'', compared to 
the X-ray ejecta parameters which suggest the remnant is plough- 
ing into material of density 0.2 — 3 cm^^ (Decourchelle et al. 2001, 
Badenes et al. 2006). The difference in the density estimated for the 
surrounding material from CO/X-ray abundances is a factor of 60 
(rather than the three orders of magnitude quoted in Tian & Leahy 
(201 1)). However, the density of the Hi cloud they detect at -43 — 
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Figure 7. Three-colour FIR image of Tycho's SNR with PACS data at 70, 
100 and 160 /xm. The area shown is 23.2' X 27.2'. 

57kms^^ is comparable to the density obtained from X-ray ob- 
servations and could therefore be tracing gas interacting with the 
remnant. Using a flat galactic rotation model (with standard con- 
stants) places this Hi cloud at d = 3.4 kpc which is consistent with 
the distance of Tycho's SNR. 

Ishihara et al. (2010) detect warm dust features in the north 
of the remnant with AKARI and suggest these are linked to the CO 
structures in the northeast. They also find enhanced dust emission 
in the northwest region of the shell and since they find no corre- 
sponding CO feature, they suggest this component could originate 
from SN dust. Without the long wavelength data, they were un- 
able to constrain the amount of cold dust arising from SN material 
within the blast wave. 

The three colour Herschel PACS image is presented in Fig.|7] 
As in Kepler, the remnant is clearly visible in blue, and the interstel- 
lar structures become clearer at the longer PACS and SPIRE wave- 
lengths (Fig. [8). Cold, dusty clouds peaking at 250 /^m are seen to 
the immediate east and west of the remnant and have temperatures 
typical of the ISM. The bright SPIRE interstellar structure to the 
immediate north of the remnant in Fig. [8] is highly suggestive of 
emission originating from gas and dust swept up by the expanding 
SNR. The structure is similar in appearance to the submm cavities 
blown out by the winds from warm OB stars (e.g. Schneider et al. 
2010) with filamentary fingers pointing inwards to the centre of the 
ionising radiation (in this case the SNR). In the higher resolution 
PACS data (Fig. |7j, we see three small, warm sources embedded 
within the SPIRE clouds (see i]4.2| for further discussion). The com- 
pact blue source to the north (Fig. [TJ is the Cepheid Variable star, 
V* AS Cas. 

4.1 The multi wavelength view 

The emission from warm dust at 70 /im in the north has two bright 
components in the NE and NW (Fig.[8]l; seen also in the AAj4/?/data 
at wavelengths out to 160 /im. This corresponds spatially to peaks 
in 24 fim emission as well as the X-ray continuum which traces the 
forward blast wave. The longer wavelength Herschel-SPW£. data 
do not spatially correlate with radio or X-ray emission from the 



remnant. In Fig. |9l we compare the PACS 70 pm emission with a 
three-colour X-ray image (0.95-1.26, 1.63-2.26 and 4.1-6.2keV). 
The emission from warm dust at 70 /im is bounded within the for- 
ward and reverse shocks of the SNR at 251" and 183" respectively 
(Warren et al. 2005) and we do not see FIR emission towards the 
interior of the remnant. The 250 /im emission (Fig.[8]l, however, ex- 
tends far outside the shock boundary and spatially correlates with 
the low-level structures and the peaks seen in the CO map (Fig.llOt. 
There doesn't appear to be any agreement between the peaks at 70 
and 250 /^m though we see emission at both FIR wavelengths where 
there is also CO emission. Comparing with the CO molecular gas, 
the interstellar clouds seen by Herschel close to the remnant in the 
north (labeled NW and NE) and to the east, match structures seen 
at —61 and —65 kms^^. In the southeast and southwest, two faint 
CO clouds are seen in the same range of velocity, both with associ- 
ated 70 ^m and 250 ^m emission. 

The correlation between the 70 ^m emission and the X-rays 
emitted at the outer edges of the SNR is striking. The hard X-ray 
emission (4.1-6.2 keV) is filamentary in structure with thin rims ex- 
tending outside of the circular edge both in the north and south-east 
side of the remnant. The peaks in the warm dust trace the peaks in 
the hard X-rays. The FIR emission does not spatially coincide with 
the soft X-rays (Fig.|9} with peaks in the soft X-ray lagging behind 
the 70 /im peaks by 27" in the northwest, though the soft emission 
at the edges of the outer blast wave do correlate with the FIR (re- 
gions A' and 'B' in Fig.|9}. The peaks in the radio emission also lag 
by ~ 16". The ejecta material does not spatially correlate with the 
IR emission, with the peak in the northwest lagging behind the FIR 
by 23". The correlation of warm dust with the hard X-ray emission 
at the edges of the remnant is slightly surprising due to the large 
fraction of non-thermal emission expected from shocks in this en- 
ergy regime. However, dust is seen where soft X-rays (from thermal 
emission due to the blast wave sweeping up interstellar material) 
and hard X-rays overlap. Lee et al. (2010) showed that the hard 
X-rays also contain a significant fraction of thermal emission, we 
therefore conclude that the warm dust is spatially correlated with 
the swept up material. The lack of 70 fj,m emission in the SW re- 
gion of the remnant, where the hard X-rays are bright, may suggest 
that the CO (Fig.[TO) and 250- /^m structure in this region orginates 
from a foreground cloud which is not associated with the remnant; 
this would mean very little interstellar dust is being swept up on 
this side (in agreement with the faster expansion rate of the SNR 
observed in this region). The hard X-rays in the west-southwest re- 
gion are produced by non-thermal synchrotron emission revealed 
by the presence of X-ray 'stripes' due to tangled magentic fields 
on small scales (Eriksen et al. 2011). Here then the hard X-rays 
are dominated by non-thermal processes and the lack of warm dust 
in this region further supports the hypothesis that we only see dust 
in the shocked, swept-up ISM where we find thermal continuum 
X-rays and Ha emission - this is highlighted in Fig. II II 

There is no correlation between the FIR and X-rays for the 
outermost ejecta material which has broken through the reverse 
shock region out to the very edge of the SNR (as labeled in Fig.|9l(. 
We conclude that there is no evidence to suggest that the warm dust 
is associated with the ejecta. 

4.2 The spectral energy distribution 

The IR-submm SED towards Tycho's SNR is shown in Fig. [T2l 
The fluxes estimated from Spitzer and Herschel within an aperture 
of radius 251" centred on the remnant, are plotted with their errors 
(Table [3j the zero-level offsets were applied to the maps before 
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Figure 8. Multiwavelength montage of Tycho's supernova remnant. From left to right: fop/ ^^CO( J = 1 — 0) data integrated from —68 < t) < — 53lcms 
18cm radio data and Spf'fter MIPS 24 fim. Middle: Hersche! PACS and bottom: SPIRE. The images are centred at RA = 6.3308°, Dec = 64.1372° (J2000.0; 
yellow cross) with size 9' X 9'. The location of the forward shock (radius 251") is indicated with the white circle. 
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Table 3. Fluxes in Jy measured towards Tycho's remnant within the 251" aperture defined in Fig. [8] References Douvion et al. 2001; '' - Ishihara et al. 
2010; AKARI fluxes measured at wavelengths: - 65/im; '* - 22.9/tm; - 90/im. ^ Measured on WISE data at 22 /im. 
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Figure 11. Zoomed in region on the eastern limb of Tycho's remnant: left: hard X-Ray middle: Ha emission and right: 70 /^m emission. The blastwave is at 
radius 251" and extends inwards to 222" (enclosed by solid red circles). The reverse shock at 183" is indicated by the red dashed line. In general, the hard X- 
ray emission is a combination of non-thermal and thermal emission from the interaction of the supernova and the surrounding material. The spatial agreement 
with the Ha emission suggests that the hard X-ray emission is dominated by thermal radiation in this region. The image is centred at RA = 6.4718° and 
Dec = 64.1372° with radius 72' X 72". 




Figure 9. Three colour X-ray image of Tycho's SNR using Chandra archive 
data showing soft thermal X-rays from ejecta and continuum (0.95 — 
1.26 keV: red) and hard X-rays tracing the both thermal and non-thermal 
continuum components (4.1 — 6.1 keV: blue). The green X-rays trace the 
ejecta emission (1.63 — 2.26 keV). White contours are 70 /im PACS emis- 
sion (7, 9, 10, 13 mjy/pixel). Two regions where the ejecta is protruding 
out of the reverse shock out to the edge of the forward shock are labeled. 
Labels A' and 'B' indicate regions where FIR emission is coincident with 
soft and hard X-rays. 



fluxes were measured - §23). As the synchrotron emission in Tycho 
has been steadily decreasing with a decay rate of 0.5 % yr~^ (Klein 
et al. 1979), the radio fluxes were corrected for this decay rate to 
represent the same epoch as the Herschel observations. The radio 
spectral slope exponent in the literature ranges from a — —0.52 
(Katz-Stone et al. 2000) to -0.61 (Green 2001); for a = -0.61 
the contribution of synchrotron to the 500 /^m flux is 4Jy (17%), 
for a — —0.52 this becomes 5 Jy. We use the former spectral index 
in this work since this describes the observed radio fluxes (Fig.ll2t. 




RA [:200CI) 

Figure 10. Comparison of the warm and cold dust with Herschel (contours) 
and molecular gas (inverted colourscale). Black contours trace 70 ^m PACS 
emission from warm dust (0.4, 5, 6, 7, 9, 10, 13, 16, 20 mJy/pixel). 
White contours trace the cold dust emitting at 250 fim 
(20, 50, 100, 150, 200, 250, 300, 400 m Jy/beam). The red solid 
circle traces the forward shock wave at radius 251", the red dashed circle 
traces the reverse shock boundary at radius 183" (Warren et al. 2005). The 
forward shock is at radius 251". Following Ishihara et al. (2010), the major 
interstellar CO clouds to the north of the remnant are labeled NW and NE. 
Note that the peaks in the 250 fim emission from cold dust are spatially 
coincident with the um'elated molecular material and not the warm dust 
in the SNR. The contribution to the 250 ^tm emission from synchrotron 
emission is < 6 per cent. 
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Figure 12. The IR-radio SED of Tycho's SNR. The previously pubhshed 
fluxes from the literature (Table |2) are shown with errors as vertical bars 
(IRAS = blue diamonds; WISE = blue triangles; AKARI = purple triangles 
from Ishihara et al. (2010)). (Note that the latter study used an aperture of 
diameter 525" compared to the 504" in this work.) The solid black curve is 
the sum of the two-temperature thermal dust components fitted to the IR and 
submm fluxes (dot-dashed and dashed lines respectively) at temperatures 
T„ = 90 K and Tc = 21 K. The red curve shows the total SED from the 
non-thermal and thermal components. The synchrotron power law fit to the 
radio data is shown by the dot-dot-dot-dashed line. 

The submm flux within the aperture defined by the outer blast wave 
is therefore well above the expected synchrotron emission. 

As described in i|3.2| we fit the IR-radio SED with two- 
temperature modified blackbodies and a power law spectrum. The 
best fit parameters from the SED are = 90K, Tc = 21 K, 
P = 1.5 and dust mass A/„ = 8.6 x 10"^ and Mc = 4.0 Mq. 
The median parameters and confidence intervals from the bootstrap 
method are listed in Table |4] The temperature maps are shown in 
Fig. [13] The warm dust is clearly seen within the remnant at the 
same location as the 24 fj,m emission and the outermost radio emis- 
sion. The cold dust map highlights a dust structure on the eastern 
edge which appears to overlap with the warm dust (Fig.[T3](b)) and 
appears to be swept up material. If this structure is at the distance of 
Tycho with Td ~ 21 K and absorption properties similar to inter- 
stellar dust, the dust mass would be ~ 2 Mq, well above the mass 
of heavy elements available to form dust. 

As with Kepler, we find no evidence for a separate cool dust 
component (at ~ 30 — 40 K) in the SNR when subtracting the 
hot dust traced by emission at 24 /^m from the 70 and 100 /im maps. 



4.3 An aside on Triggered Star Formation 

Returning to the three FIR-bright, small clumps seen in Fig.jS] one 
may jump to the conclusion that the remnant has triggered star for- 
mation in the compression of material as the blast wave moves out- 
wards. Triggered star formation has long been suggested as a viable 
star forming mechanism (e.g. Elmegreen & Lada 1977), where an 
expanding Hll region compresses gas between the ionization front 
and the shock front, causing fragmentation into gravitationally un- 
stable, star-forming cores. Supernova-triggered star formation has 
been inferred through evidence of supernova-shocked molecular 
gas (e.g. Reach & Rho 1999, Reynoso & Magnum 2001). The FIR- 
bright clumps seen by Herschel are detected at 24 /xm, suggesting 
these are star-forming cores rather than protostars. The tempera- 




(b) 



Figure 13. Dust temperature maps derived from fitting a modified greybody 
to the 24-500 /jm maps for the area surrounding Tycho's supernova remnant 
with SED components separated into (a) warm dust and (b) cold dust. The 
colourbar shows the temperature in K. The solid black circle indicates the 
forward shock radius at 251". 



tures we obtain from fitting their SEDs from 24-500 /im suggest 
each core has a warm and cold dust temperature component with 
53 — 57 K and ~ 12 K with dust masses ranging from 1.4 — 2.8 Mq 
for the two smaller cores and 14 Mq for the largest and bright- 
est core. These masses were estimated using standard interstellar k 
values and (assuming they are at the distance of Tycho) indicate the 
formation of intermediate- and high-mass stars. 

The star-forming clumps lie outside of the average shockfront 
boundary of 251" (4.6 pc) from the centre but the nearest, faintest 
clump (Fig. [8]( is consistent with the distance reached by the out- 
ermost blast wave in this region at 260" (as traced by the hard 
X-rays). Indeed, the travel-time for the shockfront to reach these 
star-forming cores at radius r ~ 5 — 6 pc from the centre is given 
by t = 2r/5u (assuming the SNR is in the Sedov phase where 
V is the expansion velocity - v ^ 4700 km s"'^ - Hayato et al. 
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2010). The shockfront would therefore be expected to reach the 
observed star-forming regions at t = 425 — 510 yrs. However, a 
simple timescale argument rules out any triggered star formation 
by Tycho's SN since the gravitational collapse of a cloud requires 
timescales of 10^ yrs, with the proceeding protostellar phase need- 
ing approximately 10^ yrs. It is possible that the stellar wind of the 
progenitor or the binary companion was responsible for creating 
a low density bubble due to its stellar wind. Such a scenario was 
suggested by Badenes et al. (2007) to explain RCW86 (confirmed 
observationally by Williams et al. 2011) and the star-forming loop 
around the SNR G54. 1-1-0.3 (Koo et al. 2008). 



5 DISCUSSION 

5.1 A theoretical model of ejecta dust formation 

Nozawa et al. (2011) have developed a theoretical model to follow 
the dust formation in a carbon deflagration event (Nomoto, Thiel- 
mann & Yokoi et al. 1984) with ejecta mass of 1.38 Mq. In this 
model, dust grains form early-on, but are almost completely de- 
stroyed over a timescale of 10® years due to the passage through the 
reverse shock. The total dust masses formed at 100-300 days range 
from 3 X 10"* — 0.2 Mq, depending on the sticking probability 
of the grains and/or whether CO or SiO molecules are formed in 
the ejecta (since this will affect the composition and mass of grains 
that condense out of the ejecta gas - Nozawa et al. 2003). Here, we 
use this model to compare the dust masses observed with Herschel 
and the mass expected to have formed and survived in the ejecta 
at approximately 400 years after the explosion event. (We refer the 
reader to the full description in Nozawa et al. (201 1) for details of 
the simulation). 

Fig. [14] (a) shows the time evolution of the total dust mass 
of supernova grains for a Type la using the dust formation calcu- 
lation with a sticking probability equal to one. We consider three 
ambient gas densities in which the ejecta is expanding into, where 
riH = 1,3 and 5 cm~"^ since this will affect the amount of de- 
struction of the dust grains we expect; the figure clearly shows 
that dust destruction is more efficient for higher ambient densi- 
ties. Using the relative ages of Kepler and Tycho at 410(440) years 
and assuming that the ejecta is expanding into an ambient gas 
density of 1 cm~^, the model predicts a surviving dust mass of 
~ 88(84) X 10""^ Mq: an order of magnitude more than the warm 
dust mass of 3.1(8.6) x 10~'^ Mq estimated from the Herschel ob- 
servations (Table |4ll. If the remnant is expanding into denser mate- 
rial, as certainly suggested for both the Kepler and Tycho remnants, 
the surviving dust mass is consequently reduced with 26(23) and 
17(14) X 10"^ Mq of dust surviving at 410(440) years for ambient 
densities of 3 and 5 cm~"^ respectively; the surviving dust mass is 
still too high to be consistent with the observations. 

The comparison between the dust mass estimated from simu- 
lations and that observed with Herschel implies that even if all the 
warm dust observed here was produced in the ejecta, the model still 
overestimates the amount of dust expected in the ejecta. One pos- 
sible explanation for this discrepancy is that the formation of dust 
grains in the ejecta is less efficient than the model e.g. the sticking 
coefficient needed to form the grains is lower, or that the radioactive 
heating from nickel and cobalt is higher, therefore inhibiting the 
formation of dust (Nozawa et al. 201 1). It is also possible that dust 
destruction via sputtering may be more efficient than the model, 
for example due to smaller grain sizes formed and/or higher am- 
bient densities; the Herschel observations would then suggest that 
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Figure 14. (a) The variation of dust mass with time since the supernova 
event in a Type la remnant based on the model of Nozawa et al. (2011). 
The dust mass with age are compared for different ambient densities of 
n = Icm"'^ (black), 3cm~'^ (blue) and 5cm~'^ (red). The vertical Hnes 
indicate the time since explosion for Kepler and Tycho. (b) The distance 
reached by ejecta dust with time since explosion for the Kepler SNR. The 
location of the forward (FS) and reverse shocks (RS) are indicated in red 
with the FS matched to the observed location at 110 — 120". The distance 
traveled by the dust grains are shown in blue (C) and black (Si). The grey 
vertical line marks 410 years and the horizontal line marks a distance of 
5.3 X 10^* cm. Note that the Si grains are destroyed just after 400 years. 

dust destruction occurs efficiently on timescales of ~400 years in 
Type la ejecta. Finally, using a lower mass absorption coefficient in 
our estimate of the observed dust masses would also alleviate this 
discrepancy. 

Similarly, we see no observational evidence for cool or cold 
dust associated with ejecta material as predicted using the dust for- 
mation model. The lack of cold ejecta dust has important conse- 
quences on our understanding of the formation of iron-rich dust 
grains. Type la ejecta are prime candidates for producing iron-rich 
dust as they synthesise an order of magnitude more ^®Ni (which 
eventually decays into Fe) than Type-II SNe. Furthermore, there 
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is a well known problem in reconciling the predicted ejecta abun- 
dances of iron using theoretical models (e.g. Thielemann, Nomoto 
& Hashimoto 1996) with observed abundances which can differ by 
more than an order of magnitude. It has long been suggested that 
one way to reconcile this issue is to 'hide' the missing iron in dust 
grains. The Nozawa et al. (2011) model predicts that any newly 
formed iron-rich dust grains (which would form in the innermost 
regions of the ejecta) would not yet be swept by the reverse shock 
at 400 days and would still be cold. Since Herschel does not de- 
tect any significant (~subsolar) mass of cool dust associated with 
the ejecta, this strongly suggests that neither Kepler nor Tycho pro- 
duced significant amounts of iron grains. These observations sug- 
gest that the formation of dust cannot explain the discrepancy be- 
tween the observed and expected iron abundances in Type la super- 
nova ejecta (see also Hamilton & Fesen 1988). 



5.2 The origin of the warm dust component 

Given that we observe less dust than expected to be formed in the 
ejecta itself, could a non-supemova source of dust explain the IR 
emission? One can estimate the mass of dust (Table g) we would 
expect to be swept-up from circumstellar and interstellar mate- 
rial by the SN blast wave using simple arguments. Kepler is not 
yet in the Sedov phase (Decourchelle & Ballet 1994; Kinugasa & 
Tsunemi 1999; Cassam-Chenai et al. 2004) which suggests that 
the remnant is still dominated by the dynamics of the ejecta and 
has swept-up a mass similar or less than that ejected in the ex- 
plosion. The interstellar density at the location of Kepler's rem- 
nant (~ 470 pc out of the galactic plane) is expected to be low 
at riH ~ 10"^ - 0.5 cm~^ (Vink 2008). Estimating the volume 
encompassed by the SN blast wave, the swept up interstellar gas 
mass would be ~ 10~^ and hence the swept up dust mass is 
< 10~* M0 (with gas-to-dust ratios of 100 — 150 - e.g. Devereux 
& Young 1990). This is two orders of magnitude lower than the 
dust mass observed in Kepler. 

As suggested by Blair et al. (2007) and Douvion et al. (2001), 
the origin for the IR emission in Kepler could be swept up circum- 
stellar dust. X-ray emission arises from this shocked dense circum- 
stellar component with mass ^ 1 M© (Kinugasa & Tsunemi 1999; 
Blair et al. 2007; Chiotelis, Schure & Vink 2011). Given standard 
gas-to-dust ratios, this mass could easily explain the warm dust 
masses seen in the IR. In this scenario, the explosion is single- 
degenerate and the progenitor has a relatively massive evolved 
companion, indeed Chiotellis, Schure & Vink (2011) suggest the 
companion was a 4 — 5 AGB star. Is the dust mass we ob- 
serve here consistent with the expected mass loss from an evolved 
companion star? Given the typical timescale for an AGB mass- 
loss phase of ~ 10® yr and typical dust mass-loss rates of Md ~ 
IQ-io - 10"* M0 yr"^ (Meixner et al., 2004; Groenewegen et al. 
2007; Lagadec et al. 2008), we would expect < 0.01 M© ejected 
during this phas^f]. This is consistent with the mass of warm IR- 
emitting dust in Kepler and supports either the explosion of an 
intermediate mass WD-AGB binary system (e.g. Chiotellis et al. 
2011) or the single progenitor 8M0 AGB-star thermonuclear ex- 
plosion suggested by Reynolds et al. (2008). 



^ Since some of the circumstellai'/interstellar dust swept up by the super- 
nova shock wave will also be destroyed via sputtering, the masses estimated 
from these simple swept-up volume arguments are therefore upper limits on 
the actual swept-up dust mass we'd expect. 





Parameters 
(K) 


Tc 
(K) 


Dust mass 
(10-3 Mq) (M0) 


Kepler 


1.5 

1.6 ±0.2 


84 


19 


3.5 

o 1 +0.8 
•^•-^-0.6 


2.2 

n 1+0.8 
^•-^-0.5 


Tycho 


1.5 ±0.2 
1.5± 


90 

90I7 


21 

2ltl 


8.6 


3.7 

t; 1+0.9 
^•-^-0.6 



Table 4. The best-fit parameters (top row) for dust towards and within Ke- 
pler and Tycho's SNRs and the median pai'ameters (bottom row) estimated 
from bootstrapping the fluxes within their en'ors. We quote the median from 
the distribution of values obtained from the bootstrap analysis with the 68% 
per cent confidence interval quoted for the error. The dust mass absorption 
coefficient, k, used to estimate the warm dust mass in Kepler's and Tycho's 
SNR is K70pm = 2 kg" ^ . For the cold dust component towards the 
SNRs, we use Ksoo^m — 0.1 m kg 



In Tycho's remnant, it is unlikely that the swept-up dust orig- 
inates from the CSM as the companion star was inferred to be a 
low mass sub-giant (Lu et al. 2011). Some of the dust could have 
been created during the prior AGB phase of the WD star; using 
the argument above, we could again expect to see 0.01 M0 of dust 
ejected in this phase which is consistent with the measured dust 
mass for Tycho. However, there is no evidence for the presence 
of a dense CSM surrounding Tycho. The ambient interstellar den- 
sity as measured from X-ray and radio evolution is often quoted as 
TiH < 0.2 cm-3 (Dwarkadas & Chevalier 1998; Cassam-Chenaf 
et al. 2004; Katsuda et al. 2010) but recent simulations suggest 
3 cm~3 (Badenes et al. 2007). We would therefore expect a swept- 
up dust mass'' of up to 0.013 — 0.2 M0 for the range of ambient 
gas densities (with gas-to-dust ratio of 120). This is consistent with 
the warm dust component we measure in Tycho's remnant. 

Recently, Tian & Leahy (2011), using Hi absorption lines, 
showed that a number of the molecular clouds thought to be in- 
teracting with Tycho are not at the location of the remnant. This 
would imply there is not enough interstellar medium surrounding 
Tycho for the remnant to sweep up the dust we observe. However, 
they identify one Hi cloud which they rule out as interacting with 
the remnant due to its density of nn ~ 10 cm-^, yet we include 
this here since this is comparable to the density of the surrounding 
material as estimated from hydrodynamical modelling (Badenes et 
al. 2006). The mass of gas in this cloud is ~ 10 M0. Furthermore, 
we identify a faint cloud overlapping with the warm dust in the NW 
of Tycho at — 55kms-^ (with gas mass ~ 1 M0 using a CO-H2 
factor of 3 X 10" cm' -2 (Kkms-^)-^) which is massive enough 
to explain the origin of the warm dust feature detected here'' and in 
Ishihara et al. (2010). 

The warm dust we observe in these remnants is therefore co- 
incident with the outermost shockfronts, indicative of swept-up cir- 
cumstellar or interstellar material. We can test this by investigating 
the location of the dust grains expected to form from the ejecta 
gas using the model of Nozawa et al. ( i]5.1| l w.r.t the forward and 
reverse shocks. In Fig. [14] (b), we show the distance the freshly- 
formed ejecta grains will have reached in 1000 years with ejecta 
ploughing into an ambient medium of density nn = 5 cm^^ (cho- 
sen so that the model forward shock location for Kepler agrees with 
the observed shock radius at ~ 110 — 120"). This figure suggests 
that silicon and carbon dust grains formed in the ejecta will en- 
counter the reverse shock after 50 years, and at 400 years (due to 
deceleration from the gas drag), will be trapped at 5.3 x 10^* cm 
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from the centre; for wh = lOcm"'^, the dust would be trapped at 
4.5 X 10^^ cm. The model thus predicts that any supernova dust 
in Kepler would lie at an angular distance of 78 — 91". The latter 
distance is consistent with the interior edge of the dust emission 
seen with Herschel (Fig. [3} but not consistent with the emission 
seen out to 120". In Tycho, the observed warm dust is confined be- 
tween a thin region at 220 — 250" from the centre (Figs. [TO] and 
II U . The model suggests that for a forward shock of ~ 250", the 
supernova dust would be located at approx. 1 x 10^^ cm, i.e. at an 
angular distance of 200", again consistent with the interior edge of 
the observed warm dust emission in Tycho. 

If the warm dust was formed in the supernova ejecta, we would 
expect it to be distributed within the forward and reverse shocks 
across the entire remnant. However, in Kepler's SNR, the warm 
dust exhibits a north-south asymmetry, and in Tycho, there is a 
north-north-east asymmetry and we therefore only observe warm 
dust emission where there is also plenty of surrounding material. 
For Kepler, this is the circumstellar material in the north which the 
remnant is ploughing into as it moves at high velocities, and in Ty- 
cho, the dust originates from regions where we see dense molecular 
clouds at the edge of the remnant. This indicates that the majority 
of the warm dust is therefore swept-up. In summary, we cannot 
absolutely rule out that a small fraction of the warm dust in both 
remnants is formed in the supernova ejecta, but we propose that, 
from the spatial distribution, the warm dust in Tycho orginates from 
swept-up interstellar dust, and in Kepler, the warm dust originates 
from swept-up circumstellar dust. 

5.3 The origin of the cold dust component 

Approximately 2 — 5 of cold dust is found in the apertures 
chosen to encompass the SNR, for Kepler and Tycho respectively 
(Table|4]l. The substantial amount rules out a swept-up CSM or ISM 
origin of dust, but could this massive cold component be freshly- 
formed dust in the supernova ejecta? 

As well as using the temperature maps (Figs[6l&|13|l, we also 
searched for an excess of submm emission in the SNR compared 
to the unrelated molecular gas. The CO and submm fluxes (using 
the 250 /im) in each pixel were compared and a best-fit relation 
determined, producing, for each pixel, an estimate of the 250 
flux that corresponds to a given CO flux. The scatter in the CO and 
250-/^m fluxes for Kepler's remnant is considerable, and when the 
submm flux associated with the CO is subtracted from the 250 
map, residual submm emission remains in the northern shell (in the 
same location as the structure seen in the cold temperature map, 
Fig|6](b)); the dust mass in this feature is ~ 0.6 for Td ~ 20 K. 
The cold dust does not spatially coincide with any of the ejecta 
tracers, and therefore we conclude that the submm emission is not 
dust formed in the supernova ejecta. 

There is no obvious CO structure at the location of the cold 
dust shell seen in Kepler, though the difference between the flux at 
250 in this region and the interstellar cloud to the east of the 
remnant is a factor of 5; the expected CO emission associated with 
this faint cold dust feature is therefore at the noise level of the CO 
map. There is submm 'excess' over the molecular gas both inside 
and outside the remnant, suggesting that these dust clouds in gen- 
eral are not always well traced by the ^^C0(J=2-1) emission; this 
method is therefore not able to distinguish between dust associated 
with the foreground dust component and the remnant. 

Towards Tycho, a large mass of cold dust exists within the 
aperture encompassing the forward blast wave. Assuming that this 
dust component originates from the swept-up interstellar material 



would imply more than 600 M0 of gas has been swept up in the 
last 440 years. This is unphysical since Tycho is still in the Sedov 
phase where the swept up mass < 10 x Mojocta. The cold dust 
on the eastern edge of the remnant as seen clearly in the temper- 
ature map (Fig. I13t appears (at first glance) to be swept-up cold 
dust, but the large mass of dust in this feature (~ 2M0) again 
rules this scenario out. There are a significant number of atomic 
and molecular clouds overlapping on the line of sight, the average 
H2 column density towards the remnant (as estimated from the CO 
map) is ~ 10^^ cm^^, with clouds of sizes 3 — 6 pc at distances of 
3.4 — 5kpc, suggesting gas masses of ~ 500 — 2000 M0. There- 
fore, given the abundance of gas towards and behind Tycho (and 
spatially bounded between the SN shock radius), it is likely that 
we are observing a small amount of cold dust swept up in the ISM 
along with a significant amount of dust emission from interstellar 
clouds along the line of sight to the remnant. 



6 CONCLUSIONS 

Herschel observations of the remnants of the core-collapse SNe Cas 
A and SN1987A show evidence for a significant mass of dust spa- 
tially coincident with SNe ejecta. Using Herschel observations of 
Kepler and the Type la SNR, Tycho, we do not find evidence for 
freshly-formed supernova dust in these historical remnants. 

We confirm a previously known warm dust component in Ke- 
pler's remnant with Td — 82lg K, somewhat cooler than the dust 
previously detected by Spitzer and ISO, and revise the dust mass 
to~ (3.llo:6) X IQ-^Mo.The warm dust is spatially coincident 
with thermal X-ray emission and optical knots and filaments on the 
outer edges of the shockfront, confirming that this component orig- 
inates in the material swept up by the primary blast-wave. Given 
the mass and the density of the surrounding ISM, the warm dust 
could not have originated from swept up interstellar material, but 
is consistent with swept-up circumstellar dust, confirming the pres- 
ence of CSM surrounding Kepler. The Herschel SPIRE images are 
dominated by emission from cold foreground insterstellar clouds 
(~ 20 K) which overlap with the small dust clumps detected by 
previous ground-based SCUBA measurements. These were origi- 
nally attributed to a large mass of dust either formed in the ejecta 
and/or swept up by the circumstellar medium. The resolution of 
Herschel does not allow us to distinguish between cold supernova 
dust or intervening interstellar material, however using dust temper- 
ature maps, and a comparison with molecular gas (traced by CO), 
we find that most of the cold dust originates from foreground, un- 
related ISM. There is no evidence that any of the cold dust within 
the aperture is spatially coincident with ejecta material. 

Similarly for Tycho's remnant, we detect warm dust at 901^ K 
with mass (S.Olj g) x 10"'' Mq. Comparing the spatial distribu- 
tion of the warm dust with thermal and non-thermal X-rays from 
the ejecta and swept up medium, and with collisionally heated Ha 
emission arising from the thin post-shock edge, we suggest the 
warm dust is swept up interstellar material. This is further sup- 
ported by the significant amount of molecular gas surrounding the 
regions where warm dust is observed and the large mass of inter- 
stellar gas swept up by the remnant. We find no evidence for cool 
supernova dust at temperatures of 20 — 30 K as detected in the core 
collapse remnants Cas A and SN1987A; the lack of freshly formed 
dust in both remnants supports the view that Kepler was the result 
of a Type la explosion. 

We compare the dust masses with the recent theoretical dust 
formation model of Nozawa et al. (2011) where approximately 
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1Q~^ Mq of dust is expected to form in Kepler and Tycho. Tlie 
Si and C grains formed are expected to tiave traversed the re- 
verse shock approximately 50 years after the explosion and are then 
trapped at 400 years at a distance consistent with the interior edge 
of the warm dust emission detected by Herschel. The observed dust 
mass of the warm component is lower than the model predicts, sug- 
gesting that dust formation in the ejecta is inhibited or destroyed 
more efficiently. We do not detect a cool dust component in the 
innermost regions of un-shocked ejecta, suggesting that Type la 
ejecta are not producing substantial masses of iron and are not re- 
sponsible for the 'missing' iron mass problem in the ejecta. The 
mass of dust in these remnants is lower than that measured in core- 
collapse SNRs, though there is a possibility that a fraction of the 
cold dust observed towards both Tycho and Kepler could be hidden 
amongst the dust components along the line of sight. Resolved ob- 
servations of younger Type la remnants with low ambient densities 
(and little foreground contamination) are needed to investigate the 
evolution of supernova dust mass further. 
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